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APPLICATION FOR 
UNITED STATES LETTERS PATENT 
SPECIFICATION 

TO ALL WHOM IT MAY CONCERN: 

Be it known that we, Kurt BRAUER, a citizen of Germany, residing at 
Kirchnerweg 3, 72076 Tubingen, Germany, and Hans-Martin HAFNER, a citizen of 
Germany, residing at Friedrich-Zun-delstrasse 30, 70619 Stuttgart, Germany, have 
invented a new and useful METHOD AND DEVICE FOR DETERMINING A 
CHARACTERISTIC VALUE RELATING TO THE BLOOD CIRCULATION OF 
MODIFIED TISSUE REGIONS IN ORGANS OF LIVING BEINGS, of which the 
following is a specification. 



Method and device for determining a characteristic value for 
the perfusion of modifi ed tissue regions in organs of living 

beings 



Related Applications 

This is a continuation of copending International Patent Appli- 
cation PCT/EP02/08855 filed on August 8, 2002 and designating 
the U.S., which has been published in German as WO 03/015626 A2 
and which claims priority of German Patent Application No. 101 
40 924 filed on August 14, 2001. 

Background of the Invention 

Field of the Invention 

The invention relates to a method for determining a character- 
istic value for the perfusion of modified tissue regions in 
organs of living beings. 



2 



The invention also relates to a device for determining a char- 
acteristic value for the perfusion of modified tissue regions 
in organs of living beings. 

Related Prior Art 

In dermatology, it is known to use various examination methods, 
for example direct-light surface microscopy, epiluminescence 
microscopy, dermatoscopy etc. in order to examine pigmented 
tissue regions. These examination methods are based on the 
discovery that the morphology of different tissue regions, 
especially the microcirculation inside pigmented tissue re- 
gions, can be characterised by characteristic values which can 
be used as evaluation criteria for this tissue region. In this 
context, it is known that the capillary vessels are pathologi- 
cally modified in malignant tumors and so-called hypervascu- 
larization is formed in the edge region of the tumorous tissue, 
with capillary blood vessels growing in the direction of the 
tumor. This phenomenon is referred to as neoangiogenesis ("in- 
grown vessels" ) . 

It is known to use laser Doppler flowmeters to examine the 
perfusion in human skin tissue (Kvernmo H.D. et al . , Microvas- 
cular Research 57, pp. 298 - 309, 1999). In this context, it is 
also known to process laser Doppler perfusion signals, which 
are recorded as a volumetric flow rate as a function of time, 
specifically by means of a so-called wavelet analysis. The term 
"wavelet analysis" is intended to mean a method in which a 
time-variant signal is three-dimensionally represented, spe- 
cifically over the axes of frequency or scaling, time and sig- 
nal amplitude. Wavelet analysis is preferably used for those 



time-variant signals in which a plurality of periodic or quasi- 
periodic processes are superimposed, possibly together with 
other stochastic processes. In such a case, wavelet analysis 
makes it possible to deliberately filter out or deliberately 
suppress particular frequency or scaling ranges in which indi- 
vidual periodic or quasi-periodic signals occur. 

The three-dimensional representation of a laser Doppler perfu- 
sion signal by means of wavelet analysis will be referred to 
below as a "vasomotion field". 

In the article by Kvernmo et al. mentioned above, the perfusion 
in human skin tissue is examined as a function of the adminis- 
tering of particular vasodilating drugs. The characteristic 
values determined in the scope of that work therefore relate 
exclusively to the effect of administering these medicaments. 
Characteristic values which relate to the perfusion behavior of 
the tissue itself, particularly in comparison between pigmented 
and nonpigmented tissue regions, are not disclosed in that 
paper. 

Summary of the Invention 

It is an object of the invention to improve a method and device 
of the type mentioned at the outset, so that it is possible to 
determine characteristic values which are indeed characteristic 
of the perfusion of modified, in particular pigmented tissue 
regions in internal organs, in particular the skin. 

In a method of the type mentioned at the outset, this object is 
achieved according to the invention by the following steps: 
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a) determining a first volumetric flow rate of the blood by 
means of laser flowmetry at a first measurement point 
within the modified tissue region as a first measurement 
profile of a volumetric flow rate signal as a function of 
time ; 

b) determining a second volumetric flow rate of the blood by 
means of laser flowmetry at a second measurement point 
within an unmodified tissue region as a second measurement 
profile of a volumetric flow rate signal as a function of 
time; 

c) carrying out a wavelet analysis of each of the measurement 
profiles as a three-dimensional representation of the sig- 
nal intensity over a frequency or scaling axis and time; 

d) determining the profile of a vasomotion energy over the 
frequency or scaling axis for each of the measurement pro- 
files subjected to the wavelet analysis , the vasomotion 
energy being the integral of the signal intensity with re- 
spect to time for a particular frequency or scaling value; 

e) subtracting the profile of the vasomotion energy of the 
second measurement profile from the profile of the vasomo- 
tion energy of the first measurement profile, so as to 
form a first difference profile; 

f) integrating the difference profile in order to obtain the 
characteristic value. 



In a preferred variant of the method mentioned at the outset, 
the object of the invention is furthermore achieved by the 
following steps: 

a) determining a first volumetric flow rate of the blood by 
means of laser flowmetry at a first measurement point 
within the modified tissue region as a first measurement 
profile of a volumetric flow rate signal as a function of 
time ; 

b) determining a second volumetric flow rate of the blood by 
means of laser flowmetry at a second measurement point 
within an unmodified tissue region as a second measurement 
profile of a volumetric flow rate signal as a function of 
time; 

c) determining a third volumetric flow rate of the blood by 
means of laser flowmetry at a third measurement point at 
the transition between the modified tissue region and the 
surrounding, unmodified tissue region as a third measure- 
ment profile of a volumetric flow rate signal as a func- 
tion of time; 

d) carrying out a wavelet analysis of each of the measurement 
profiles as a three-dimensional representation of the sig- 
nal intensity over a frequency or scaling axis and time; 

e) determining the profile of a vasomotion energy over the 
frequency or scaling axis for each of the measurement pro- 
files subjected to the wavelet analysis, the vasomotion 
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energy being the integral of the signal intensity with re- 
spect to time for a particular frequency or scaling value; 

f) subtracting the profile of the vasomotion energy of the 
third measurement profile from the profile of the vasomo- 
tion energy of the second measurement profile , so as to 
form a first difference profile; 

g) subtracting the profile of the vasomotion energy of the 
first measurement profile from the profile of the vasomo- 
tion energy of the second measurement profile, so as to 
form a second difference profile; 

h) forming an average value of the first and second differ- 
ence profiles; 

j) integrating the average value in order to obtain the 
characteristic value. 

In a device of the type mentioned at the outset, the object of 
the invention is achieved by the following features: 

A device for determining a characteristic value for the perfu- 
sion of modified tissue regions in organs of living beings, 
having: 

a) a laser flowmeter for determining 

a first volumetric flow rate of the blood at a first 
measurement point within the modified tissue region 
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as a first measurement profile of a volumetric flow 
rate signal as a function of time; 

- a second volumetric flow rate of the blood at a second 
measurement point within an unmodified tissue region 
as a second measurement profile of a volumetric flow 
rate signal as a function of time; 

b) means for carrying out a wavelet analysis of each of the 
measurement profiles as a three-dimensional representation 
of the signal intensity over a scaling axis and time; 

c) means for determining the profile of a vasomotion energy 
over the frequency or scaling axis for each of the meas- 
urement profiles subjected to the wavelet analysis, the 
vasomotion energy being the integral of the signal inten- 
sity with respect to time for a particular frequency or 
scaling value; 

d) means for subtracting the profile of the vasomotion energy 
of the first measurement profile from the profile of the 
vasomotion energy of the second measurement profile, so as 
to form a difference profile; 

e) means for integrating the difference profile in order to 
obtain the characteristic value. 

In a preferred variant of the device mentioned at the outset, 
the object of the invention is furthermore achieved by the 
following features : 
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A device for determining a characteristic value for the perfu- 
sion of modified tissue regions in organs of living beings, 
having: 

a) a laser flowmeter for determining 

- a first volumetric flow rate of the blood at a first 

measurement point within the modified tissue region 
as a first measurement profile of a volumetric flow 
rate signal as a function of time; 

- a second volumetric flow rate of the blood at a second 

measurement point within an unmodified tissue region 
as a second measurement profile of a volumetric flow 
rate signal as a function of time; 

- a third volumetric flow rate of the blood at a third 

measurement point at the transition between the modi- 
fied tissue region and the surrounding, unmodified 
tissue region as a third measurement profile of a 
volumetric flow rate signal as a function of time; 

b) means for carrying out a wavelet analysis of each of the 
measurement profiles as a three-dimensional representation 
of the signal intensity over a frequency or scaling axis 
and time; 

c) means for determining the profile of a vasomotion energy 
over the frequency or scaling axis for each of the meas- 
urement profiles subjected to the wavelet analysis, the 
vasomotion energy being the integral of the signal inten- 
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sity with respect to time for a particular frequency or 
scaling value; 

d) means for subtracting the profile of the vasomotion energy 
of the third measurement profile from the profile of the 
vasomotion energy of the second measurement profile, so as 
to form a first difference profile; 

e) means for subtracting the profile of the vasomotion energy 
of the first measurement profile from the profile of the 
vasomotion energy of the second measurement profile, so as 
to form a second difference profile; 

f) means for forming an average value of the first and second 
difference profiles; 

g) means for integrating the average value in order to obtain 
the characteristic value. 

The object underlying the invention is fully achieved in this 
way . 

In fact, it has been found that the features listed above make 
it possible to determine characteristic values which permit 
virtually 100% determination of the properties of the modified 
tissue region being examined. This reliability has also been 
confirmed in a clinical study. 

According to a further object, the characteristic value may be 
compared with a reference value. 
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This measure has the advantage that it is possible to assign 
the characteristic value which is determined to clinically 
determined limit values, so as to draw conclusions therefrom in 
order to describe the state of the modified tissue region. 

Other advantages can be found in the description and the ap- 
pended drawing. 

It is to be understood that the features mentioned above, and 
those which will be mentioned below, may be used not only in 
the respectively indicated combination but also in other combi- 
nations or individually, without departing from the scope of 
the present invention. 

Brief Description of the Drawings 

An embodiment of the invention is represented in the drawing 
and will be explained in more detail in the description below. 

Fig. 1 shows a highly schematized representation of a skin 
region, in which there is a pigmented tissue region, 
while specifying measurement points selected according 
to the invention ; 

Figs 2A 

and 2B show laser Doppler perfusion signals, as can be deter- 
mined at the measurement points according to Fig. 1 as 
volumetric flow rates of the blood as a function of 
time; 
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Figs 3A 

to 3D show four vasomotion fields, as are obtained by wavelet 
analysis from laser Doppler perfusion signals; 

Fig. 4 shows the vasomotion field according to Fig. 3D, on an 
enlarged scale and in order to explain further details; 

Fig. 5 shows the graphical profile of the vasomotion energy 
over a scaling axis for a malignant pigmented tissue 
region, in comparison with the surrounding tissue; 

Fig. 6 shows a representation similar to Fig. 5, but for a 
benign pigmented tissue region; 

Fig. 7 shows a bar chart to explain a clinical study carried 
out according to the invention. 

Description of a Preferred Embodiment 

In Fig. 1, 10 denotes a skin region of a human in a highly 
schematized way. Within the skin region 10, there is a modified 
tissue region 12, which is a pigmented tissue region in the 
exemplary case being represented. The pigmented tissue region 
12 may be benign or malignant. Benign manifestations include 
freckles, liver spots and the like (melanocytic nevi ) , while 
malignant manifestations include the various types of skin 
cancer (melanomas). The skin region 10 stands for an example of 
any organ of a living being, particularly an internal organ. 

In Fig. 1, 14 denotes a first measurement point at the centre 
of the pigmented tissue region 12, 16 denotes a second measure- 
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ment point at the edge of the pigmented tissue region 12, and 
18 denotes a third measurement point in the surrounding, non- 
pigmented tissue region. The second measurement point 16 is not 
absolutely necessary for carrying out the method, although it 
is evaluated in preferred variants of the method. 

At the two measurement points 14, 18, or the three measurement 
points 14, 16 and 18 , measurements are then carried out by 
means of a laser Doppler flowmeter. Such laser Doppler flow- 
meters are known to the person skilled in the art, and are sold 
for example under the device reference "Periflux System 5000" 
by Perimed, Stockholm, Sweden. 

In laser Doppler flowmeters, a laser beam is guided into the 
vessels to be examined. Owing to the Doppler effect, the light 
reflected by the moving blood corpuscles has its frequency, 
i.e. color, shifted as well as its intensity. In this way, it 
is possible to carry out measurements of the volumetric flow 
rate of the blood. Measurement profiles of a volumetric flow 
rate signal as a function of time are then obtained. 

Figs 2A and 2B represent two such measurement profiles 20 and 
22 for the measurement points 14 and 18 in Fig. 1. It can be 
seen clearly that a larger signal amplitude is obtained at the 
first measurement point 14 and a smaller signal amplitude is 
obtained at the third measurement point 18. The measurement 
Profiles 2 0 and 22 also reveal that a plurality of periodic and 
quasi-periodic processes seem to be superimposed on a stochas- 
tic process. The periodic or quasi-periodic processes are in 
this case the vital physiological functions of the human body 
(respiration, pulse etc.). These processes change in a known 
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way with time, so that the measurement profiles according to 
Figs 2A and 2B per se are not very informative. It can, how- 
ever, be inferred therefrom that the pigmented tissue region 12 
seems to be perfused substantially more than the surrounding 
nonpigmented tissue region. 

In order to make the recorded measurement profiles 20 and 22 
more informative, a special wavelet analysis is carried out 
according to the invention. As already mentioned above, wavelet 
analysis is a tool known per se in signal processing to convert 
a time profile into a three-dimensional representation with the 
axes of time, frequency or scaling and signal intensity, for 
measurement profiles of the type available here. From the wave- 
let analysis representation, the various periodic, quasi- 
periodic and stochastic processes can then be identified better 
and, if need be, filtered out or suppressed. 

Figs 3A to 3D present four such representations, the abscissa 
being the time axis t and the ordinate being the frequency or 
scaling axis s. The third axis of signal intensity is produced 
in the representations of Figs 3A to 3D by differential darken- 
ing (or coloring). 

The representation in Fig. 3A corresponds to a reference pat- 
tern for benign tissue, the representation in Fig. 3B corre- 
sponds to a central measurement for benign tissue, while Figs 
3C and 3D represent the corresponding reference and central 
measurements for malignant tissue. 

For further explanation, the representation of Fig. 3D is shown 
again on an enlarged scale in Fig. 4. The right-hand half of 



14 



Fig. 4 represents four typical regions p 0 , p lf p 3 and p 5 . The 
region p 0 , which is not in fact shown in the actual image , 
corresponds to the respiration (slow frequency), which merges 
through a transition region p x into the region p 3 of the heart- 
beat (medium frequency), which is followed further above after 
another transition region by the region p 5 of the reflected 
pulse waves (high frequency). 

This provides a characteristic pattern with periodic bands, 
which are assigned to the respective periodic processes. 

In order to interpret this wavelet analysis, so-called scaling 
levels are defined on the frequency or scaling axis s, one of 
which is indicated in Fig. 4 by 30 in the region p 3 of the 
heartbeat. By progressing along the scaling level 30 in the 
abscissa direction through the diagram according to Fig. 4 and 
integrating the (preferably squared) signal values with respect 
to the time axis, then a value is obtained for each scaling 
level 30 which will be referred to below as a "vasomotion en- 
ergy E 11 • The vasomotion energy E can be plotted for the various 
scaling levels s, as is done in Figs 5 and 6. 

In Fig. 5, a first profile denoted 32 and a second profile 
denoted 34 of the vasomotion energy E are plotted against the 
scaling level (or frequency axis) s. The first profile 32 be- 
longs to a malignant melanoma (for example corresponding to the 
wavelet analysis according to Fig. 3D, i.e. according to Fig. 
4), while the second profile 34 corresponds to the surrounding 
nonpigmented skin region. It can be seen clearly that there is 
a significant difference between the profiles 32 and 34. 
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For comparison, Fig. 6 shows corresponding profiles which were 
obtained from the examination of likewise pigmented but benign 
tissue regions in comparison with the surrounding nonpigmented 
tissue. A third profile 38 corresponds to the nonpigmented 
surrounding region in this case, while the fourth profile 36 
corresponds to a melanocytic nevus. Apart from the region be- 
tween s = 1 and s = 2, the profiles are not so significantly 
different as in the case of Fig. 5. 

In order to explain how a reliable characteristic value can be 
determined from such measurements, reference will again be made 
to the representations according to Figs 2A and 2B in conjunc- 
tion with Fig. 1. 

According to a first variant of the invention, the profile of 
the vasomotion energy E corresponding to the second measurement 
profile 22 (surrounding nonpigmented tissue) is now subtracted 
from the profile of the vasomotion energy E of the first meas- 
urement profile 24 (centre of the pigmented region 12). 

The difference profile obtained in this way is then integrated. 
This integral is the intended characteristic value, which will 
be referred to below as AE. 

According to a second variant of the invention, however, the 
profile of the vasomotion energy E corresponding to the second 
measurement profile 22 (edge region) is subtracted from the 
profile of the vasomotion energy E of the third measurement 
profile 24 (surrounding nonpigmented tissue), and the profile 
of the vasomotion . energy E of the first measurement profile 20 
(centre of the pigmented region 12) is likewise subtracted from 
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the profile of the vasomotion energy E, also of the third meas- 
urement profile 24. 

The average value is then formed from the difference profiles 
obtained in this way, and this average value is integrated. 
This integral is the intended characteristic value, which will 
be referred to below as AE. 

Fig. 7 represents the result of a clinical study, which was 
carried out on 61 patients. Measurements corresponding to the 
representations in Figs 1 and 2 were carried out on these pa- 
tients, and these measurements were then evaluated in the way 
described above. When determining the characteristic value AE, 
they were normalized to a range of between 0 and 1 for better 
clarity. 

After the measurements, samples were taken from the patients in 
the respective pigmented tissue regions 12, and these were 
histologically examined. 

In Fig. 7, the bars 40 show the results of those patients for 
whom benign manifestations of the pigmented tissue region were 
found in the course of the histological examination. For exam- 
ple, about 30% of these patients had a characteristic value AE 
of between 0.3 and 0.4. From the bars 40, it was possible to 
determine a first distribution curve 42 for these patients 
having benign pigmented tissue regions. 

The bars 44 represented as dark in Fig. 7, however, belong to 
those patients for whom malignant pigmented tissue regions were 
established subsequently in the histological examination. From 
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the bars 44, it was possible to derive a second distribution 
curve 4 6 for these patients. 

As can be seen clearly from Fig. 7, the two distribution curves 
42 and 46 are separated relatively distinctly from each other 
on the AE axis. It is therefore possible to specify a limit 
value 50 which, with the selected normalization to values of 
between 0 and 1, lies at AE = 0.75. 

In the scope of the present invention, the wavelet analysis is 
preferably carried out by using the following concepts: 

The measurement profiles 2 0 and 2 2 according to Figs 2A and 2B 
can be analytically described as follows: 

J samp 

where £(t n ) represents the measurement points at the times t n , 
and a constant sampling rate f samp was selected. The vasomotion 
field V at the time t T on a scaling level s is then obtained 
from the wavelet representation of the time profile as: 



nej 1..NJ 
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The wavelets <j> s are derived from a so-called mother wavelet <t>. 
For analysis of the vasomotion data in the scope of the present 
invention, the following mother wavelet was found to be par- 
ticularly suitable: 



0(t) = cos (17.3932 • 0 e~' 2 .. 

This mother wavelet describes a local oscillation over about 
ten periods. The wavelets are derived from the mother wavelet 4 
in the following way: 



</> s {t)=4^<p{st). 



Vs normalizes the wavelets <|> s independently of the scaling 
level s . 



The scaling levels are subdivided into so-called "octaves" and 
"voices " : 



s -> s a = 2 alNiv ° K " ) mit cr e {o..N (voices) M ioctaves) } . 

N (OC taV 03) ±s number Qf oct - aves in VaSOm0tion field/ 

NfmcM) is tne number of voices per octave. This leads to the 
so-called continuous wavelet transformation. 
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The voices and octaves of the scaling levels are closely re- 
lated to the frequencies of the Fourier spectrum. For better 
understanding, the scaling values are indicated in units which 
are as close as possible to the Fourier frequencies. Five dif- 
ferent scaling ranges are then obtained, as already explained 
above with reference to Fig. 4: 

PO: s CT =0.15 - 0.31 (respiration) ; 

PI: s a = 0.31 - 0.62 (transition from respiration to heart- 
beat); 

P2: s a = 0.62 - 1.20 (transition from respiration to heart- 
beat) ; 

P3: S* » 1.20 - 2.50 (heartbeat); 

P4: s a = 2.50 - 5.00 (transition from heartbeat to reflected 
pulse waves ) ; 

P5: s a = 5.00 - 10.00 (reflected pulse waves). 
Overall, the vasomotion field V is then given as: 



ne{\..N} 
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Such vasomotion fields V are represented as two-dimensional 
printouts with gray-value- or color-coded amplitudes in Figs 3A 
to 3D and 4 . 

For further evaluation of the signals/ the vasomotion field 
energy E S(J is determined for each scaling level s a , by summation 
over the time t. 

e =Tv 2 

ztgN 



For comparison of the different vasomotion fields V, in par- 
ticular between healthy and pigmented tissue, the energy dif- 
ference or characteristic value AE is determined as follows: 

AE=C * I ^(healthy) _g( pigmented) I 

€ 4 | s a s <j y 

S,CT 

C= Z I E ( s h J althy) +E { s p ^ mented) |=> AE <=[0,1 J 



